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Abstract. We present a model for nebular emission in star forming galaxies, which takes into account the effects of dust 
reprocessing. The nebular emissions (continuum emission, 54 H and He recombination lines and 60 nebular lines from UV to 
IR) have been computed with CLOUDY and then included into GRASIL, our spectrophotometric code specifically developed 
for dusty galaxies. The interface between nebular emission and population synthesis is based on a set of pre-computed Hn 
region emission models covering a wide range of physical quantities (metallicity, density, geometry and number of H i. He i and 
11 ionizing photons). These quantities are fully adequate to describe the emission properties of the majority of star-forming and 
starburst galaxies. Concerning the extinction properties of normal star forming galaxies, we are able to interpret the observed 
lack of correlation between the attenuation measured at Hor and in the UV band as a consequence of age selective extinction. We 
also find that, for these galaxies with modest SFR, the ratio FIR/UV provides the best constraints on the UV attenuation. The 
accurate treatment of lines and continuum in dusty galaxies also allows to deal with different SFR estimators in a consistent way, 
from the UV to radio wavelengths, and to discuss the uncertainties arising from the different physical conditions encountered 
in star forming galaxies. We provide our best estimates of SFR/luminosity calibrations, together with their expected range of 
variation. It results that SFR derived through Ha, even when corrected for extinction using the Balmer decrement, is affected 
by important uncertainties due to age selective extinction. Another remarkable result is that SFR from UV luminosity corrected 
by means of the ratio FIR/UV has a small uncertainty. Finally, our model provides a calibration of SFR from radio luminosity; 
its value differs from estimates from other works, but we are able to reproduce the observed FIR/radio ratio. These results are 
relevant to estimates of the contribution of disk galaxies to the cosmic SFR at z < 1 . 
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J> 1. Introduction wealth of information is coming from infrared spectral region. 

r Both star formation and a possible nuclear activity contribute 

^r' Massive stars in young stellar generations leave their signature . ^ • u . -.u j-a . • • 

rN , J b b = to infrared luminosity, but with different emission properties 

*H as UV, ionizing photons and supernova (SN) explosions. Thus, i . i i aaa t-. u- e j-a j- 

rrt ' \ ' f (Rigopoulou et al. 19991 The combination of different indi- 

■ - - 1 ongoing star formation in galaxies is traced by the reprocess- . / u i j j • f j • • v 

° ° ° J r cators (such as optical, near- and mid-intrared emission lines, 

ing of stellar radiation (line and continuum emission from the j ■ e jnAou j u t^u t e 

° ^ mid-infrared PAH bands, shape of the continuum from the far- 
surrounding ionized gas, thermal emission by dust grains) and • 1 -.u A- 1 »U N U 1J»U 11 » Ut ■ 

° ° & / ultraviolet to the radio wavelengths) should then allow to obtain 

by synchrotron emission of SN accelerated electrons. The hy- ^ • • . ^ . f j ui az-^m 

•'a fair picture of the star formation process and possible AGN 

drogen recombination lines have been widely used as tracers of » u * j i * i • 

° ^ contribution for a wide class of active galaxies, 
the current star formation rate (SFR), because of the direct pro- 
portionality between their intensity and the number of living Two models have been recently published in which spec- 
massive stars. In the absence of dust, nebular emission models ^r^l synthesis and ph oto-ionization codes ai-e coupled. Chariot 
predict an almost constant ratio of Ha to H/3 for a large variety & Longhetti (2001 henceforth CLOl) combined a popula- 
of environments (e.g. Osterbi-ock .l989» : so, the dust absorp- tion synthesis code with the photo-ionization code CLOUDY 
tion can be estimated by comparing the observed ratio to the (Ferlandt2001 1, including an approximate prescription to es- 
theoretical one (Balmer decrement method). tim^te the dust absorption. According to CLOl, the presence 
However, the discovery by IRAS of powerful and highly °f ^ust increases the uncertainty of standard SFR estimators 
dust enshrouded starbursts (Soifer et al. 1986), namely the ^^^^^ as Ha or [On] luminosities, from a factor of a few to 
Luminous and Ultra Luminous Infrared Galaxies, has under- s^^ei'^l decades. Only the simultaneous consideration of sev- 
lined the necessity to model nebular emission by consistendy ^ral other lines and spectral features can reduce the uncertainty 
taking into account dust reprocessing. Moreover, an impressive wi*in a factor of a few. In a similar way, Moy et al. (|200T| 
hereafter MRFOl) interfaced the evolutionary synthesis model 
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these models the dust processes are partially treated (extinction 
but not emission) using a simple screen approximation. ' 

However, growing evidence has been collected in the re- 
cent past showing that a physical understanding of dust effects 
in galaxies requires the inclusion of different environments, ar- 
ranged in a rather complex geometry. A sophisticated treatment 
tends to become more and more important as the obscuration 
gets higher 

In this paper, we propose a new method to compute nebu- 
lar emission in star forming galaxies. The method was imple- 
mented in the spectrophotometric synthesis model GRASIL^ 

The main advantage with respect to previous treatments is 
that GRASIL already provides a sound treatment of all the as- 
pects of dust reprocessing, since it makes use of a geometry 
which is much closer to reality than a screen between the stars 
and the observer. Dust modeling is thus related to real phys- 
ical parameters, describing the different distributions of dust 
and stars in molecular complexes and diffuse components, and 
their age dependence. 

One of its important features is that it takes into account 
that, since stars are born in dense environments (the molecular 
clouds) and progressively become less obscured, the relative 
geometrical arrangement of dust and stars depends on the age 
of the stellar generation considered. Granato et al. (2000) have 
shown that this age-dependent extinction can explain the dif- 
ferences between the observed attenuation laws in normal and 
starburst galaxies (see also Poggianti et al. 2 (TOT). 

GRASIL has been shown to reproduce the UV to radio con- 
tinuum SEDs of galaxies, at low and high redshift, and in dif- 
ferent evolutionary stages (Silva et al. 1 19981 Franceschini et al. 
[TW8. Granato et al. .2000^ Rodighiero et al. EDOOi Granato et 
al-IMD- 

The paper is organized as follows. In Sect. |2 we describe 
how we modeled stellar radiation and the computation of neb- 
ular emission from the single Hii regions. Sect.|3lexplains the 
calculation of nebular emission from galaxies and in the fol- 
lowing sections our model is applied to the case of normal star- 
forming galaxies. In particular in Sect.l^the different methods 
to estimate the attenuation are discussed. In Sect. [S] different 
SFR estimators, together with their uncertainties, are presented. 
In Sect.|6lwe discuss the use of IR nebular lines. Thus, in Sect. 
0the main results of the paper are discussed, and then summa- 
rized in the last section. Finally, some technical details of Sect. 
12.2. ll andBl are described in appendixIXlandlHI 

2. Population synthesis with gas and dust 

In this section we summarize the main features of our popu- 
lation synthesis code GRASIL. More details can be found in 
Silva et al. (|T998 ) and Granato et al. ( ^0001 . GRASIL repre- 
sents galaxies by means of two main components characterized 

' Additional details on these models are given and compared to our 
method in Sect.|3| 

^ The GRASIL code and updated information on it can be found at 
the URL: 

jhttp : //web . pd . a stro. it/granato/grasil/grasil. h tml | 
(Silva et al. I1998t , and then it was used to study the properties of 
normal star-forming galaxies. 



by different geometries: a spheroidal component (the bulge) 
and a disky component. Dust, that may be present in both com- 
ponents, is divided in two phases: i) dense molecular clouds 
(MCs), where star formation is active, and ii) diffuse medium 
(or cirrus). 

Young stars are assumed to be born into MCs, and to leave 
them progressively as their age increases. As a consequence, 
the fraction of light of young simple stellar populations (SSPs) 
radiated inside MCs is a decreasing function of SSPs age, pa- 
rameterized by the "escape time". Thus, the light of young stars 
will be attenuated by both the MCs and the cirrus, while older 
populations will only be affected by dust in the cirrus compo- 
nent. The molecular cloud is modeled as a thick spherical shell 
of dense gas (and dust) around a central point source, represent- 
ing all the stellar content of the cloud. The time dependence 
of the escape fraction gives rise to the age-selective extinction 
because younger stellar generations are more attenuated than 
older ones. 

Massive stars are supposed to ionize the surrounding 
medium and to give rise to the corresponding line and contin- 
uum nebular emissions. This radiation comes generally from 
within MCs and it is accordingly extinguished. However, it is 
worth noticing that, when the escape time is short enough, a 
significant number of ionizing photons can arise from star gen- 
erations outside MCs and, consequently, we will also consider 
Hii regions extinguished only by the cirrus component. To 
compute nebular emission we have used the photo-ionization 
code CLOUDY (version 94), as detailed in the following sec- 
tions. 

Since our interest is focused on the star formation pro- 
cess and the obscuration of star-forming regions, we do not 
include in our model the emission due to an Active Galactic 
Nucleus (AGN). We neglect the ionization due to UV radia- 
tion from post-AGB stars, since in presence of even a modest 
star formation activity their contribution is very low (Binette 
et al. ll994> . Moreover, we neglect the contribution of shocks 
produced by SN explosions, that is typically low (Kewley 
et al. .2001 ) . Finally, we do not include line emission from 
photo-dissociation regions and diffuse warm neutral/low ion- 
ized medium, although some fine structure IR lines are effi- 
ciently produced in these media. As a consequence their lumi- 
nosities in our model will be only lower limits; this problem 
will be discussed in Sect. 12.31 

The final output of our model is a complete and detailed 
spectrum of star-forming galaxies, from far-UV to the radio 
wavelengths (an extension to X-ray band is in progress, Silva 
et al. in preparation), including stellar absorption features, neb- 
ular emission, dust and PAH emission. 

2.1. Stellar Radiation 

SEDs of stellar generations (SSP, simple stellar populations) 
have been computed by following the prescriptions outlined in 
Bressan et al. (1994) and in Silva et al. (1998 l. These SEDs 
cover a wide range in age and metal content. They can be com- 
puted for an arbitrary initial mass function (IMF) and allow the 
use of different atmosphere models, from the low-resolution 
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(but covering a wide parameter space) Kurucz-Lejeune mod- 
els (Kurucz 11993 1 Lejeune et al. I1998I I. to the intermediate- 
resolution models of Pickles (1998 1 and Jacoby et al. ( 1984 1. In 
the latter two cases, which are derived from observed stars, the 
corresponding fluxes have been extended into the unobserved 
region by means of the Lejeune et al. ( 1998 1 models. The use 
of higher resolution models is particularly useful when dealing 
with emission lines superimposed to absorption features of the 
intermediate age populations. As for the most massive stars, we 
have adopted the atmospheric models by Schaerer et al. ( 1996 ) 
for mass-losing blue supergiants, and the models by Schmutz, 
Leitherer & Gruenwald 119921 in the Wolf Rayet (WR) phase. 
It will be stressed below that the emission from H ii regions 
has been computed in a way that is almost independent of the 
detailed shape of the adopted SEDs of stellar populations. 

2.2. H II regions model 

The emission spectrum from a single Hii region depends on 
two main ingredients: the SED of the ionizing star cluster and 
the properties of the excited gas. For a given stellar evolution 
scenario, the SED is determined by the IMF and the total mass, 
age and metallicity of the cluster, as well as by the adopted 
model atmospheres. On the other hand, the important proper- 
ties of the excited gas are the density, the chemical composition 
and the geometry. In a galaxy the situation is more complicated 
and star forming regions possibly have diflferent ages and differ- 
ent metallicities. In starbursts, star clusters of diff'erent ages and 
metallicities may coexist within the same star forming region. 
To compute the line emission intensities, one should consider 
the spectrum of the ionizing source, provided by the recent star 
formation history, and use a photoionization code with suitable 
values of the gas parameters. To reduce the computing cost, one 
can use a pre-built library of line intensities, corresponding to 
different ages of the ionizing SSP and different metal contents 
of both the SSP and the gas. This still requires different libraries 
for different masses of the ionizing SSPs, assumed IMFs and 
model atmospheres. It becomes particularly time consuming in 
applications requiring a large number of models. 

Our approach has been to pick out the physical parameters 
which actually affect the emission properties of H n regions. 
The analysis described in the following section shows that the 
emission line spectrum of an H ii region with fixed gas proper- 
ties (metallicity, density and geometry) is described with rea- 
sonable precision by only three quantities: the number of ion- 
izing photons for Hi, Hei and On (Qh, Que, and Qo, defined 
in eq.EJ- This method allows us to get rid of the particular SSP 
model and IMF. In fact, different ionizing sources that provide 
the same values of Qu, Qhs, and Qo, will produce the same 
emission line spectra, within a reasonable accuracy. 

Thus, we computed a library of photoionization models as 
a function of Qu, Qh&, and Qq. When estimating the actual 
line emission due to a given stellar population, we compute Q^, 
Qiif., and Qo from the corresponding SED, and then interpolate 
the value from the above library. 

The H II region library can be freely retrieved through the 
GRASIL web page (see Sect.[lJ- 



In the following, we will describe the adopted procedure 
and test its accuracy. 

2.2.1 . Ionizing spectra 

In order to have photoionization models as a function of Qn, 
Que, and Qo, we approximate real SEDs with piece-wise 
blackbodies, in a way such that the space of Q parameters is 
covered as much as possible. 

The Lyman continuum spectra of typical young stellar pop- 
ulations show two prominent discontinuities, at 504.1 A and 
227.8 A, corresponding respectively to the ionizing energies 
for He I and He ii. At wavelengths smaller than 227.8 A there 
are very few photons. Between the Lyman break, Hei break, 
and He ii break, spectra can be reasonably well represented by 
blackbodies. Thus, when computing the libraries, we used the 
following functional representation of the ionizing SEDs: 



r for ^ < 227.8 A 

AueByiTHe) for227.8A</l<504.1 A 

AiiBy{TH) for 504.1 A < ^ < 911.76 A 

AniBv(7'ni) for911.76 A < /I 



(1) 



where By(T) is the Planck function at temperature T. This func- 
tion depends on the 6 quantities A^e, T^e, Au, T^, A„i and T^i- 
However, we found that the SEDs of young stellar populations 
can be well approximated by writing these 6 quantities in terms 
of only three parameters. These are the numbers of ionizing 
photons for H i. He i and O ii (Qu, Qh&, and Qo)'- 



Qh = T^dv , gHe = T^dv , go = 

Jvv «^Vh.- •^Vn 



hv 



-dv. 



(2) 



where vh, vhe and vq are the photoionization threshold frequen- 
cies for respectively H i. He i, and O ii (the last one correspond 
to Ao - 350.7 A). The relationship between the Q values and 
the parameters of analytical spectra is described in appendixIXl 
We built a library for a grid of values of Qu, Q^f., and Qo, 
and for different assumptions on the gas density, metallicity and 
filling factor (e). When the filling factor is different from 1, the 
gas in the H ii region is supposed to be divided in small clumps, 
and e is defined as the ratio between the volume occupied by 
the clumps and the total volume of the Hii region. The gas 
density refers to the hydrogen density inside the clumps. Hii 
regions are assumed to be spherical and ionization bounded, 
with a constant density along the radius, and with a covering 
factor of 1 . The abundance of elements respect to H are from 



McGaugh ( II991> . These abundances are relative to gas phase, 
thus they already account for the fraction stored in dust (deple- 
tion). 

The characteristics of dust inside H n regions and its effects 
on the emerging emission are still poorly known. However, it 
is reasonable that the main effect is that a fraction /j (typically 
~ 30%, see DeGioia-Eastwood. I1992> of ionizing photons is 
absorbed by dust. To approximate the consequence, one can 
multiply the nebular emission of a dust-free model by a factor 
1 - /d. Thus, absorption by internal dust has not been consid- 
ered in H II region models presented here. Moreover, including 
dust inside photo-ionization models would introduce new free 
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parameters, and the library could be coherently used only in 
some specific cases. 

It is worth noticing that the emission properties of the ion- 
ized gas (i.e. the ratio of different line intensities) at a fixed den- 
sity do not depend directly on the geometry of the gas, while 
it shows a dependence on the ionization parameter U, that ex- 
presses the number of ionizing photons per particle. For a thin 
spherical shell at distance R from a central ionizing source, the 
ionization parameter is [/ = gH/(4^^^c«H); it can be shown 
(see appendixlEl that the volume averaged ionization parame- 
ter {U) for a sphere, with density hh and filling factor e, is 



3af(H,r)/3eH«He 



{U) 



Ac 



2\l/3 



47r 



(3) 



where a^(i{, T) is the recombination coefficient of hydrogen 
in the case B. Thus, in our geometry (U) is proportional to 
(Qunne^y^^ while for the plane parallel geometry U oc Qu/nu. 

As a consequence, our results remain also valid for different 
geometries and/or filling factors, provided that we renormaUze 
Qii to maintain the same value of (U). 

The computed lines are 54 H and He i recombination lines 
and 60 nebular line of other elements with wavelength from UV 
to FIR. The complete list of lines and the space of parameters 
covered by the library can be found in the web page. 

2.2.2. Nebular continuum emission 

The nebular continuum emission can be easily evaluated us- 
ing the Stromgren theory for the case B approximation (see 
Osterbrock 1989j and added to the emerging spectrum. 
We get 



yv(H,r) rv(He,r) 



aB(H, T) 



■aB(He, T) 



(4) 



where yv(He, 7") and yv(H, T) are the total emission coeffi- 
cients and aB(He, T) and q'b(H, T) are the recombination coef- 
ficients, respectively of helium and hydrogen. In the former co- 
efficients we take into account recombination, free-free emis- 
sion of H n and He ii and the Lya two photons decay emission. 
Free-free emission is computed as described in Bressan et al. 



J2002> . while other coefficients are from Burgess & Summers 
( fT976l . Aller (1984i, Osterbrock ( fT989t . Ferland ( 1 980 1 and 
Nussbaumer & Schmutz (I1984> . The Heiii contribution to the 
nebular emission is very small and can be neglected. 

2.2.3. Accuracy 

Before proceeding further, we need to analyze the errors intro- 
duced by adopting the analytic SEDs (equation^ instead of 
the real SSP with the same Q values. 

A first problem is that our analytical spectra neglect pho- 
tons with energy higher than the Hen ionization threshold 
(hereafter we will refer to these photons as high energy pho- 
tons). This does not allow a realistic evaluation of the lumi- 
nosity of the He ii lines. However, nebular He ii lines are su- 
perimposed to the generally more intense photospheric lines, 
emitted in the stellar atmospheres of WR stars (Conti, .1991. 
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Fig. 1. Differences between Hii region models excited by SSP 
spectra with and without high energy photons. Crosses are for 
difference smaller than 1%, empty squares for difference be- 
tween 1% and 5%; circled symbols refer to very low luminos- 
ity lines. Gas and star metallicity is 0.008; the stellar mass is 
3 ■ lO-* Mo. 



Schaerer et al. ll999> . We also cannot compute line intensities 
of highly ionized elements, such as O iv or Ne iv. These lines 
are typically either not observed, or very faint in pure star form- 
ing galaxies, so that they are often interpreted as a signature of 
AGN activity. The impact of neglecting high energy photons on 
other lines is quantified by comparing the calculated emission 
of Hii regions excited by SSP spectra with and without high 
energy photons, normalized in order to have the same Qh- The 
results (for some ages and metallicities) are shown in Fig.^for 
a Salpeter IMF between 0.15 and 120 Mq. The difference in 
the two cases is always < 10%, and typically < 1%; the highest 
difference appears at SSP ages at which the WR stars produce 
their maximum flux in high energy photons. We conclude that 
the lack of high energy photons is not very important, at least 
for the lines considered here. 

The comparison between line intensities of Hii region 
models excited by an SSP spectrum and those excited by a 
parametric spectrum with the same values of Qu, Que, and Qo 
is shown in Fig.|2 Typical differences in the emission lines are 
lower than 10%, with the larger ones arising at lower values 
of Qii, i.e. older ages. Most of them are due to discrepancies 
between the SSP and parametric SEDs at wavelengths longer 
than the Lyman break. Indeed, photons with wavelengths corre- 
sponding to Lyman serie lines can be efficiently absorbed and 
re-emitted in optical lines. When the ionizing flux decreases 
with respect to the non-ionizing one, this effect can be impor- 
tant. Notice however that, when computing the lines emitted by 
a combination of H ii regions excited by clusters with different 
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Fig. 2. Differences between Hii region models excited by SSP 
spectra and parametric spectra with the same Qs. Empty trian- 
gles represent differences between 5% and 10%, filled triangles 
represent differences greater than 10%. 



ages, the most important contribution arises from the youngest 
populations. 

As noticed in previous works (see e.g. Stasinska et al. 120011 



or Rubin et al. l2001t . the photo-ionization models cannot be ac- 
curate when the gas metallicity exceeds the solar value, for sev- 
eral reasons. As the metallicity increases, fine structure infrared 
lines of metals (that depend little on temperature) dominate the 
cooling processes, but their transition probabilities have not yet 
been well determined. On the other hand, optical collisionally 
excited lines, that depend strongly on the electronic tempera- 
ture, become uncertain because the temperature is regulated by 
infrared lines. Furthermore, the thermal instability of the gas 
affects the line intensity when the temperature is low and the 
metallicity is high. Finally, a further source of uncertainty at 
high metallicity is due to the poorly known depletion of metals 
into dust grains. In order to minimize this problem, we used 



a set of abundances measured by McGaugh J1991> in the gas 
phase. 

In order to check the flexibility of our procedure, we com- 
pare in Fig. 13 the results obtained by repeating the above test 
when adopting a very different IMF, namely a top-heavy mod- 
ified Kennicutt's IMF (d log A^/d log m = -0.4 for 0.15 Mq < 
m < 1 Mo and d log A^/d log m = -1 for 1 Mq < m < 120 Mq, 
Kennicutt. ll983> . The differences between Hii region models 
excited by SSP spectra and analytical ones with the same Q val- 
ues are not larger than in the case of the Salpeter IMF, despite 
the very different slope of the IMF. 

In conclusion, our analytical spectra can reproduce with a 
good accuracy the nebular emission properties of H ii regions 
excited by star clusters with different IMF, age and metallicity 



e [%] 

Fig. 3. As Fig. |3 but with SSPs with a top-heavy modified 
Kennicutt's IMF. 



by only making use of the three quantities Qh, Que, and Qo- 
Obviously this accuracy could be improved by adding some 
other quantity to describe the SED with more details, but this 
would dramatically increase the number of H n region models 
used for the library. 

2.3. C and O infrared lines 

Ci, C II and O i fine structure infrared lines are produced not 
only in H ii regions, but also in warm neutral and ionized in- 
terstellar medium and photodissociation regions (PDR). [Oi] 
63.2yum and 145.5//m, [Ci]369yum and 610^(m infrared lines 
have upper energy levels respectively of 228 K, 326 K, 62.5 
K and 23.6 K (Kaufman et al.'1999"(, so they can easily be pro- 
duced in neutral medium. Because carbon has a ionization po- 
tential ( 11 .26 e V) which is lower than H, the C ii ion is present 
in PDR and in neutral medium illuminated by far-UV stellar 
radiation. Indeed, the [Cn]157.7yum line is the most important 
coolant of warm neutral medium. The relative contribution of 
different media to these lines is still a matter of debate (Heiles 
1994 Malhotra et al. 2001 1. Since our model does not include 
PDR or neutral gas emission, the luminosities predicted for 
these lines must be taken as a lower limit. 

2.4. Comparison witii observed H ii galaxies 

In Fig. 0] we compare our Hii library with observations of a 
sample of H n galaxies (crosses) from Dessauges-Zavadsky et 
al. ( 20Wl, a revision of the Terlevich's catalogue (Terlevich et 
al. 1991 ). Hii galaxies are typically defined as galaxies with 
spectra very similar to those of Hii regions; these objects are 
interpreted as bursting dwarf galaxies. In the figure we show the 
dependence of models on (U) (or Qu via the equation |3j and 
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Fig. 4. Diagnostic diagram Log([Oiii]/HjS) vs. Log([S ii]/Hq'). Crosses are observed Hii galaxies from Dessauges-Zavadsky et 
al. (2000 1. The three solid closed lines show the regions occupied by Hii region models from our library with three values of 
Qne/Qn (0.5, 0.275 and 0.05). The arrows show shifts at increasing {U}, Qhc/Qh, or Qo/Que- Log((f/)) ranges from -0.867 
to -3.534, and Qo/Que from 0.13 to 0.63. The models have a gas metallicity of 0.004, while the closed dotted line shows the 
position of models with solar metallicity and Qhb/Gh - 0.275. All the models have a hydrogen density of 10 cm"^. The thick 
long dashed line separates the region occupied by H ii galaxies (on the left) from the AGN region, as empirically defined by 
Veilleux & Osterbrock (il987i( : the thick dot-dashed line delineates the same separation computed by Kewley et al. (1200 U . 



on the hardness of the ionizing SED, expressed by Qne/Qn and 

Go/eHe. 

The library covers the region of the [Oiii]/Hj0 vs. [S iil/Ha 
diagnostic diagram populated by observed H ii galaxies. Some 
models fall in the region of the diagnostic diagrams occupied 
by AGNs; these models are produced by the hardest spectra in 
the library (Ghc/Gh - 0.5), harder than the spectra that can be 
buih with ti-ue SSPs. 



3. Nebular emission from a star-forming galaxy 

When considering a star-forming galaxy as a whole, we need 
to integrate the emission lines resulting from the effects of dif- 
ferent stellar populations. This cannot simply be obtained by 
calculating the ionizing photon fluxes Qh, Qhs, and Qo of the 
integrated spectrum and by getting the corresponding line in- 
tensities from our interpolation tables. Indeed, the intensity of 
metal Unes also depends on the hardness of the ionizing spectra 
and on the ionization parameter To obtain a more realistic de- 
scription of the nebular emission, it is then necessary to model 
the formation and evolution of the population of H ii regions, 
each one characterized by a different intensity, hardness and 
ionizing parameter This may be accomplished by splitting the 
recent SF history in subsequent episodes of suitable duration, 
and by computing their separate contributions to the emission 
lines. 



The total emission in the line I (Ei) at the epoch of obser- 
vation r, may be written as 



(5) 



where Nunj is the number of Hii regions that have formed in 
the time interval [T-tj+i , T-tj], q j are the corresponding 
ionizing photon fluxes, Zgas is the current metallicity of the gas 
and El is the emission from the single H ii region. 

Nhuj is obtained by assuming that each Hii region is illu- 
minated by a single cluster of total mass M*, so that 



J''*' ^'(T - t)dt 



M" 



(6) 



where *P is the SFR. Q*^ j for the single H ii region at different 
ages can be obtained by 



Xr Xr - t)(SAt,Z(T - t))lhv)AtAv 



(7) 



where we have emphasized the dependence on the metallicity 
of the ionizing spectra. Similar expressions hold for ., and 

The emission properties of the j population depend on the 
ratios QJjg^/QJj^ and Q*qj/QIi^j, that are constrained by the 
age of the j population, and on the average ionization parameter 
of the H n regions {U}*. 
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The correspondence between the integrated Lyman contin- 
uum of the population j (expressed by A^hh.j XQhP and (f/)* is 
determined by the stellar mass M* of ionizing clusters through 
equations |6] and 



1 2/3 



An 



-2\ 



1/3 



(8) 



As shown in Fig.|4] Hii galaxies exhibit a big variation of the 
ionization parameter value, which means that the value of stel- 
lar mass M* and/or of the filling factor are subject to substantial 
variations. Observations (e.g. Kennicutt 1984') suggest that M* 
varies approximately from 1000 to 10* M© and e from 0.1 to 
0.001. 

The duration of the subsequent episodes in which we split 
the SF history should be chosen by considering the evolution 
of the shape of the ionizing spectra. A lower limit is set by 
considering that each H ii region will be illuminated by all the 
stars formed within the finite formation time of a typical star 
cluster, that we assume to be 1 Myr (see e.g. Fuente et al. l2001> . 

As remarked in Sect. 12.2.11 the presence of dust inside Hii 
regions can be important. Thus one can multiply the nebular 
emission obtained from eq.|5]by a factor 1 - fa- 

Finally, the emission lines are extinguished in the same way 
as the stellar populations that produced them (see Sect.|5J. 

In galaxies with normal dust content, almost all the pho- 
tons in the resonant Lyes' line are reprocessed via two-photon 
decay and dust absorption and re-emission. In our model we 
assume that all Lya photons decade via two-photons emission, 
and then these photons are possibly absorbed by dust. Anyhow, 
their contribution to the IR luminosity amounts to < 10%. Only 
in the case of very low dust content, the relative contribution 
of absorbed Lyes' photons increases respect to the contribution 
from absorbed stellar continuum and may be significant. A de- 
tailed treatment of the Lya transfer should also include the ef- 
fects of winds and outflows, which strongly reduce the resonant 
scattering. 

We conclude this section by comparing our method with 
other approaches found in literature. 

Fioc & Rocca-Volmerange ( 19971 computed Qh from the 
integrated spectrum and then analytically obtain the HjS line in- 
tensity. Other hydrogen and metals lines are then rescaled in 
fixed proportions. Moy et al. ( 2001 1 have improved the previ- 
ous model by coupling the integrated spectrum with CLOUDY. 

A more detailed analysis of emission properties in star 
forming galaxies is by CLOl. These authors obtained the emis- 
sion from the galaxy by summing the line intensities of sep- 
arate SSPs of given constant average metallicity, weighted by 
the corresponding SFR. CLOl also have a more detailed model 
for dust absorption (see Chariot & Fall l2000l that accounts for 
an age selective extinction. 

Compared to MRFOl, our model accounts for the coexis- 
tence of different populations of H ii regions, while at variance 
with CLOl, it quickly computes the line emissions for different 
choices of M* and gas density. Also, the H ii regions library can 
be used for different IMFs, and we relax the hypothesis that the 
metallicities of the excited gas and of the ionizing stellar popu- 
lation are the same (but a non-negligible difference is expected 



only in extreme situations). Finally, our library can be easily 
used for other population synthesis models. 

4. Attenuation properties of normal star forming 
galaxies 

In star forming galaxies, the intrinsic UV luminosity is directly 
related to the star formation rate; on the other hand, the UV 
flux is heavily affected by extinction. Therefore it is important 
to understand how dust absorption affects the UV luminosity 
and, as a consequence, the SFR estimates. 

In this section we discuss several methods to estimate the 
attenuation in normal star forming galaxies by considering the 
UV, optical emission lines and FIR properties. We will show 
that observations require the extinction of different stellar pop- 
ulations to vary with age, with younger populations suffering a 
larger extinction. 

Notice that we will use the term attenuation when it refers 
to the amount of light lost at a give wavelength from a complex 
system (i.e. a galaxy) and the term extinction for the light lost 
from a simple source (as a star) with dust along the line of sight. 

4.1. Estimation of UV attenuation 

Attenuation in local star forming galaxies can be derived in 
several ways. 

A common method is to relate the UV attenuation to the 
attenuation in Ha. Calzetti (>1997,) suggested that the attenua- 
tion of the stellar continuum is only a fraction (~ 0.44) of the 
attenuation of the ionized gas. By extending this assumption to 
the UV while adopting a suitable extinction law, it is straight- 
forward to obtain (at A^y = 2000 A, e.g. Buat et al. l2002> 



12000 



1.6A 



Ho- 



(9) 



The gas attenuation at Ha is derived from the Balmer decre- 
ment: 



Alia = 1-086- 



1 



- 1 



\jna Liy} 



(10) 



where ej^a - Tuj^/rna is derived from the adopted extinction 
law and jn/j/ jua is the ratio of Hy6 and Ha emission coefficients 
(typically assumed 2.87 at an electronic temperature of lO'* K, 
see Osterbrock 1989 1. 

An alternative method to derive the UV attenuation is to 
consider the reprocessing of star light into the infrared emis- 
sion. Meurer et al. ( 1999 ) (see also Calzetti et al. 2000i as- 
sumed that the energy re-emitted by dust in a galaxy is provided 
primarily by the UV flux of young star populations. They re- 
lated (eq. 10 of their paper) the ratio between the FIR and UV 
fluxes (now A = 1600 A), FpiR/^'ieoo, to the UV attenuation 
by: 



Aieoo = 2.5 log 



0.84Fi6oo 



+ 1 



(11) 



where the infrared flux FpiR is defined as the flux in the [40- 
120] jum interval derived from the 60 and 100 jum IRAS bands 
(Helou et al.'1988'l, and Fieoo - Aimi ■ /leoo (in W cm"^). Note 
that this relation was derived for starburst galaxies. 
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In the same framework, Buat et al. <1999> proposed a re- 
lation between the UV attenuation and the Ffir/^uv ratio 
(/luv - 2000A), which is suited for normal star forming galax- 
ies: 



Table 1. Values of the relevant parameters used to model nor- 
mal star forming galaxies. Ref. model: see Sect.|5l 



i2000 



= 0.466 + log 



1+0.433 


log 


/ ^FIR \ 


\F20m) 







(12) 



Buat et al. J2002I hereafter B02) compared the Auv derived 
from the ratio Ffir/Fuv (eq. ll2> with the attenuation suffered 
by Har in a sample (called SFG sample) of normal star form- 
ing galaxies. The SFG sample consists of 47 spiral and irregu- 
lar galaxies in nearby clusters. They were observed in the UV 
(^uv = 2000 A) w ith the SCAP, FOCA and FAUST instru- 
ments (Boselli et al.'2001V in the optical (Gavazzi et al. 2002 ), 
and in the FIR by IRAS. The galaxies were selected to have 
EW(Ha)> 6 A in order to minimize the errors in Ylct and 
fluxes. Galaxies with Seyfert activity were excluded from the 
sample. Metallicity of the galaxies in the sample ranges from 
~ Z0/4 to ~ 2Z0. 
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Fig. 5. On the ordinate, A2000 derived from the ratio Ffir//^2ooo 
(using eq. I12> : on the abscissa, Ahc, from Balmer decrement. 
Filled circles show SFG sample data. Lines connect models 
(see Sect. 14. 2> with the same Mg and star formation history, 
but different tmc- Solid lines: Mq - 10" Mq, Vsch - 0.3 Gyr"', 
Tinf = 18 Gyr, solar metaUicity; dashed lines: Mg - 10'° Mq, 
Vsch = 0.7 Gyr"', Tint = 6 Gyr, solar metallicity. To empha- 
size the dependence on metallicity, we show the first model 
(solid lines) for Tesc = 3 Myr but Z - 0.04 (dot-dashed line) 
or Z = 0.004 (dotted line). For simplicity only models with 
= 45° are shown here. 



In Fig.|5lwe report the attenuation in Wa versus the value of 
Auv derived from the ratio Ffir/Fuv for the SFG sample, as in 
Fig. 2 of their paper. As pointed out by Buat and collaborators, 
the two quantities show a lack of correlation, contrary to what 
is expected from eq.|9l 



Param. 


All models 


Ref. model 




6 -18 Gyr 


12 Gyr 


Vsch 


0.3 - 0.7 Gyr-' 


0.3 Gyr-' 


Mg 


10'"- 10" Mo 


10'" Mo 


Z 


0.004 - 0.04 


0.02 


Tesc 


1 - 9 Myr 


3 Myr 


Tmc 


0.1-1.25 


0.5 
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0° - 90° 
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0.3 


0.3 



4.2. Auv simulated normal star-forming galaxies 

In order to give an interpretation of this observational problem, 
we simulated a set of disk galaxies by exploring the space of 
parameters appropriate for normal star forming galaxies. We 
computed the star formation history, gas fraction and metal en- 
richment with our chemical evolution code (Silva et al. ll998> . 
In order to point out the dependence of the results on metallic- 
ity, we considered different values for the metallicity of stars 
and gas obtained from the chemical code. 

The parameters that regulate the star formation history in 
our models are the baryonic mass of the galaxy (Mg), the gas 
infall time scale (Tinf), and the star formation efficiency Vsch of 
the assumed linear Schmidt law. The age of the galaxies has 
been set to 10 Gyr The parameters that regulate the attenua- 
tion are the escape time Tesc, the optical thickness of MC at 1 
jum (tmc), and the orientation of the disk galaxy with respect 
to the celestial plane (9 - 0° means face-on models). The dust 
(extinction and emission) properties are similar to the galactic 
one (see Silva et al. '1998' for more details). The dust/gas ra- 
tio is assumed to be proportional to the metallicity. Emission 
lines are computed for different gas densities rtu, filling factors 
e and stellar masses of clusters M*; however we concentrate 
on H recombination lines that do not depend on hh, M* or e. 
Table [2 summarizes the values of the parameters used in our 
computations. 

The models are compared with data in Fig.|5l the attenua- 
tion at Ha is derived from the Balmer decrement (eq.^| where 
we assume = 1 .47), while the attenuation in UV has been 
derived by using eq. For simplicity, we only represent the 
most dusty (solid lines) and less dusty (dashed lines) models; 
all other cases range between the two. Models refer to a 45° 
inclination; face-on models show a slightly lower attenuation, 
while edge-on models have larger attenuations and scatter. 

The models cover quite well the location of the observed 
galaxies in this diagram and confirm a real lack of correlation 
between the UV attenuation, as derived from the Ffir/^uv ra- 
tio and AHff. 

In order to clarify the origin of this scatter and to identify 
a good estimator of UV attenuation, we contrast the intrinsic 
UV attenuation - directly extracted from our models (Auv - 
-2.5 log(Luv/iuvo)) - respectively with A^a from the Balmer 
decrement and /^fir/^'uv (see Figs.|6land0. 
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Fig. 6. UV attenuation in the models versus the attenuation at 
Ha derived from the Balmer decrement. Almost vertical lines 
connect models with the same optical thickness of MC (tmc) 
and escape time (rgsc), but different SF histories and dust con- 
tent in cirrus. Only the less and most dusty models are associed 
to symbols which refer to the escape time. The thick solid line 
refers to equation|9] 



4.2.1 . Attenuation from Balmer decrement 

The scatter in Fig.|6lmust be entirely ascribed to the interplay 
between the different stellar lifetimes associated to the emis- 
sion properties and the geometry set by the critical escape time. 
In fact. Ho- is mainly produced by ionizing massive stars with 
a lifetime around 3 Myr, while UV is also emitted by less mas- 
sive and longer-living stars. We may devise the following typi- 
cal cases. 

Escape time shorter than the typical lifetime of an ionizing 
star, crosses. Independently of tmc. the models tend to define 
a relation which seems however steeper than eq.|5] This is due 
to the fact that both Ha and UV emissions are mainly produced 
outside MCs, so that they do not respond to differences in tmc 
and the attenuation is mainly due to the diffuse medium. 

Escape time longer than the typical lifetime of an ionizing 
star, triangles. Emission lines are produced only inside MCs, 
so A-na essentially measures tmc- When tmc > 0.5 the UV flux 
produced inside MCs is completely reprocessed into the IR. 
Thus, Auv saturates while Ahb still increases. This results in an 
almost horizontal displacement in Fig.|6l 

Escape time comparable to the typical lifetime of an ioniz- 
ing star, open squares. At increasing tmc, the Ha to H/3 ratio 
increases from the emitted value to a maximum value fixed by 
the H/3 emitted outside the MCs and by the sum of the Ha still 
coming from within the MCs (the attenuation is lower at Ha 
than at HfS) and from outside. Then the ratio decreases again to 
the asymptotic value fixed by the attenuation of the diffuse gas. 
This causes a behavior that is intermediate between the former 
two cases and, in particular, gives rise to the turnover shown 
at high Tmc- As a consequence, the attenuation in Ha derived 
from the Balmer decrement is underestimated. 
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Fig. 7. UV attenuation vs the ratio ^^fir/^uv; lines refer to eq. 
[TTI eq.^land our fits. 



The variation in dust content of the galaxy and/or of the in- 
clination, at constant tmc and Tesc, shifts the models in Fig.|6l 
along a constant direction, somewhat steeper than that of eq.|5J 
furthermore, a variation in metallicity corresponds to a varia- 
tion of the amount of dust which varies the attenuation due to 
the diffuse component. These effects add further dispersion to 
the data. 

Our models suggest that galaxies in Fig. [S] with Ah,, > 1.5 
mag may be characterized by escape times larger than the typ- 
ical lifetime of the ionizing stars. Models with Aua > 2.5 have 
an equivalent width EW(Ha) lower than 6A; there are sev- 
eral possible explanations for the higher observed Aua with 
large EW(Ha), such as a small burst (increasing the equiva- 
lent width), or an underestimate of H/S or [N ii] lines (giving 
higher Ah^). 

The results of models and observations described above 
point out that age selective extinction is present in normal star- 
forming galaxies and complicates the picture of extinction in 
galaxies. The simple hypothesis of screen extinction does not 
work in this kind of galaxies; different ways to compute atten- 
uation give inconsistent results because they rise from specific 
populations (ionizing and non-ionizing UV emitters) that have 
different lifetimes and live at different optical depths. 

4.2.2. Attenuation from tfie continuum 

A more robust estimation of the UV attenuation can be ob- 
tained by considering the energy reprocessed by dust in the IR. 

In Fig. we plot the intrinsic Auv against the ratio 
Ffir/Fuy. The solid line represents eq. Eland the short dashed 
line refers to eq.ll IP. Models show that, contrary to the case of 
Balmer decrement, the UV attenuation correlates quite tightly 
with the Ffir/Fuv ratio. The origin of this coiTelation is that 
the FIR flux is essentially provided by the UV emission of 



^ To convert attenuation at 1600 A of eq. 1111 into Auv at 2000 
A we assume Fieoo = ^2000 and A2000 = 0.9 • A1600 = 0.9 • 
2.5 log(FFiR/0.84F20oo + 1), as in B02. 
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young stellar populations as assumed by Meurer et al. (|1999> 
and Buat et al. ^399). Thus, the Ffir/Fuv measures the ra- 
tio between the reprocessed energy flux and the residual en- 
ergy flux in the UV. Outside the MCs, a contribution to the FIR 
comes from the optical emission of old stellar populations; this 
causes a dispersion around the average relation. 

Similar results were obtained by Gordon et al. (J2000I); their 
models show that the Auv-^^fir/^uv relationship is valid for 
different assumptions on dust properties. The dispersion found 
by Gordon et al. (.2000j is much smaller than in our simulations, 
but in their models all the stellar populations suffer the same 
extinction (i.e. no dependence on age is introduced). 

Models also show that edge-on systems tend to follow a 
different relation from what is found for face-on systems; this 
is due to the greater contribution of the diffuse medium to at- 
tenuation in the direction of the plane of the galaxy. 

Our models suggest the following relation for face-on sys- 
tems: 



A2000 = 2.5 log 



+ 1 



1 .25/^2000 

while for edge-on systems it is preferable to use 



A2000 = 2.03 log 



0.5^21 



+ 1 



000 



(13) 



(14) 



In Fig. we also compared the relations proposed by 
Meurer et al. (IT999t (eq.[TT) and by Buat et al. ii^99\ (eq.[l3 

with models. The former seems to better agree with our model, 
and lies between face-on and edge-on models, while the sec- 
ond tends to underestimate the UV attenuation (but only by 0.8 
mag in the least favorable case). 

Meurer et al. ( 1999.) showed that the observed locus in 
the diagram Ffir/Fuv vs the UV spectral index yS in starburst 
galaxies could be well reproduced by a screen dust model with 
increasing optical thickness. They obtained the following rela- 
tion between Auv and /3: 

Ai6,)o=4.43-Hl.99/J. (15) 

Recently, Bell J2002t compared the spectral index /3 with 
Auv(FIR/UV) for a sample of normal star-forming galaxies 
and found that /? is typically greater (redder) than the expected 
value for starbursts. 

We also found a loose y6-Auv correlation for disk galaxies. 
Indeed, it is worth noticing that in these galaxies the contribu- 
tion from stars outside star forming regions is quite substantial 
and can significantly affect even the UV slope. 

5. SFR estimators 

In the previous section we justified the tangled relation between 
the attenuation properties of star forming galaxies as seen at 
different wavelengths, by means of a model that accounts for 
the complex interplay between geometry, obscuration time and 
stellar lifetimes. 

Using these models, we have obtained calibrations of the 
SFR for a wide range of different observables from the UV to 
the radio regime. Only models with Tesc ^ 9 My are used here 



Table 2. Calibrations of SFR (SFR/Luminosity) from dust 
emission, in 10"^^ Mgyr^'W"' (or a: 10"^^ Moyr"'W"'A), and 
from UV luminosity (lower panel, in 10"^^ Moyr 'W 'Hz). 



xsanQ 








^ QQ 


Q 1 T 

J. LZ— / .jy 


IR 


4.98 




FIR 


8.82 


6.12-16.4 


PAH 7.7//m 


46.7 


31.3-81.5 a 


ISO LW3 


79.9 


70.9-135. 


MIPS 24 


48.0 


34.6-114. 


MIPS 70 


9.23 


6.19-17.0 


MIPS 160 


7.66 


6.15-13.7 


UV 


284.5 


103.4-622.8 


UV^ 


125.9 


92.7-184.0 


UVfir 


104.5 


90.7-114.5 



because the majority of the data in Fig.lsjcan be explained with 
short escape times, ^ 3 Myr, in agreement with other evidences 
coming from the analysis of the UV SEDs and from the number 
counts in HR diagrams of massive stars (Silva et al.2998i(. In 
order to give a useful reference value of the calibration, we 
selected a reference model whose parameters are summarized 
in Table[n This model is represented by a star in Fig.|5] Unless 
otherwise specified, all calibrations refer to a Salpeter IMF with 
minf = 0.15 M0 and mup = 120 M© and solar metallicity. 

5.1. SFR from the continuum 

Infrared luminosity is one of the most common SFR estimators. 
The interstellar dust is able to convert in IR emission a substan- 
tial fraction of the UV light emitted by young stars. Thus, the 
IR luminosity is proportional to the SFR. 

The relation between the SFR and the infrared luminosity 
is summarized in Tabled The first column indicates the ob- 
served quantity, the second column provides the value of the 
calibration for the reference model described above, and the 
last column provides the range of variation of the calibration 
among the set of models. 

The row labeled IRg-iooo refers to the total dust emission 
from 8 to 1000 fj.m, while IR refers to the infrared emission 
estimated wifli the four IRAS bands (Sanders & Mirabel[T996) 
and FIR refers to the far-infrared emission estimated from the 
60 and 100 yum IRAS bands (as described in Helou et al. 1988}. 

In the following rows of Table|2l we show the calibration of 
the specific luminosity at the peak of the PAH emission feature 
at 7.7;um (PAH 1.1 ^im), the flux in the ISO band LW3, and the 
SIRTF experiment MIPS at 24//m, 70//m and 160yum. 

The infrared estimators exhibit larger variations in normal 
star-forming galaxies than in starburst galaxies. Inspection of 
the above table shows that the calibrations can vary by up to 
80% around the reference model. In the case of normal galax- 
ies a considerable part of radiation from young stars is not ab- 
sorbed by MCs clouds, so that variations in the duration of the 
obscuration by MCs and in the optical depth of MCs produce 
an important change in IR emission. Furthermore, the absorp- 
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Table 3. Calibrations of SFR in UV band in the dust-free case, 
of metallicity. 



and for radio emission (in 10 Moyr 'W 'Hz) as a function 



z 


UV 


1.49 GHz 


8.44 GHz 


0.0008 


78.9 


59.5 


166.1 


0.004 


88.7 


65.1 


212.6 


0.008 


94.8 


66.6 


230.9 


0.015 


99.7 


67.7 


246.2 


0.02 


103.4 


68.6 


258.4 


0.03 


107.7 


69.8 


269.1 


0.04 


112.3 


71.1 


280.7 


0.05 


117.4 


72.5 


293.5 



tion of light from old populations by diffuse dust may be a sig- 
nificant source of IR radiation, thus weakening the correlation 
with the SFR. This is quite different from more powerful star- 
bursts where most of the energy produced by young stars is 
converted into IR emission, so that variations in absorption do 
not produce any important variation of the IR flux. 

The ratio between SFR and UV luminosity (/luv=2000 A) 
is presented in the lower panel of Table |2l First row refers to 
UV uncoiTected for dust extinction, the second to the UV lu- 
minosity corrected for extinction through the UV slope (Auv 
vs yS, eg. 1 151 converted to 2000 A) and the last row to UV lu- 
minosity corrected for extinction by adopting the relation Auv 
vs logCFpiR/Fuv) with eq.^land^J These calibrations refer 
to face-on models; calibrations from edge-on models have a 
much larger range of variation, with the important exception of 
the FIR-corrected estimator. 

It is important to notice that the FIR-corrected UV estima- 
tor is quite robust; in fact, as shown in sec. 14.21 the attenuation 
in the UV is tightly related to the ratio fpiR/Z^uv, for a signif- 
icant variation of the optical depth and escape time (Fig. 0. 
Consequently, by adopting the ratio Ffir/^'uv to estimate Auv, 
one may get the intrinsic UV flux and obtain a reliable SFR in- 
dicator. In summary, UV and FIR fluxes alone are fragile SFR 
indicators, while their combination is a very good estimator of 
the UV attenuation. Notice also that, because of the slope of the 
relation Auv vs /^fir/^uv, this method is not much affected by 
uncertainties in the FIR flux due to old populations. 

We also provide the dust free calibration of UV for different 
metallicities in Table |3l 

Finally, we report in Table |3l the new radio calibrations at 
1.49 GHz and 8.44 GHz, according to Bressan et al. 12002 1. 
as a function of the stellar metallicity. The emission at these 
frequencies in our model is provided by ionized gas (thermal 
emission) and SN remnants (non-thermal emission). The vari- 
ation in luminosity (at constant SFR) is due to the variation 
of the number of ionizing photons of SSPs with metallicity, of 
the temperature of the ionized gas for free-free component, and 
to the variation of SN rate for non-thermal contributions. It is 
worth noticing that the variation of 1 .49 GHz calibration with 
metalhcity is +15%. 



Table 4. Calibrations of SFR for emission lines (in 
lO-^^Moyr-'W-i, or a: in lO'^^ Moyr-'W"'. No dust: the in- 
trinsic calibrations in the case without dust. Ref. model: cal- 
ibrations for the reference model. Max value: the maximum 
value found in our models. The values between brackets in the 
forth column refer to the edge-on case. 



Line 


No dust 


Ref. model 


Max value 




0.187 


0.758 


7.94 (25.8) 


Hao 


6.08 


6.62 


15.12 (29.3) a 


Pa/3 


0.99 


1.415 


2.52 (4.37) 


Pa^o 


0.99 


0.99 


1.06(1.18) 


Paa 


0.436 


0.505 


0.676 (0.938) 


Bry 


5.69 


6.31 


7.86(10.2) 




3.28 


3.50 


4.12(4.96) 


BlQ- 


1.69 


1.721 


1.91 (2.10) 


HUQ- 


12.9 


13.12 


14.7 (16.2) 


[On] 3727 


0.127 


0.844 


25.2 (104.8) 



Observations indicate that FIR and radio emissions are 
strongly correlated over a wide range of IR luminosities, from 
star-forming to starburst galaxies (Sanders & MirabelJ996): 



q = log 



Ffir/(3.75 X lO'^Hz) 
^'i.49GHz/(W m-^Hz-') 



2.35 + 0.2. 



(16) 



From calibrations in Tableland in Table|3l for the FIR and the 
1.49 GHz, we obtain q - 2.32 for the reference model, with a 
range of variation q - 2.05 - 2.48. 

Remarkably, the 0.2 dex observed scatter is fully accounted 
in our models simply by the variation of the SFR-FIR relation- 
ship (see also Perez-Olea & Colina 1995). This suggests that 
the relationship between radio emission and star formation is 
extremely tight and implies that the underlying physical mech- 
anism responsible for it is quite homogeneous. 

Although IR and radio calibrations depend on /d (here as- 
sumed to be 0.3), their values change by less than 5% (25% for 
8.44 GHz) if /d = is assumed. 

5.2. Optical near- and mid-infrared emission lines 

H recombination lines are the only ones almost proportional 
to the Lyman continuum flux and, as a consequence, to the 
SFR; other emission lines also depend on other quantities so 
that their use as SFR estimators is quite dangerous. 

In Table 0] we report our calibrations for some emis- 
sion lines. The second column shows the intrinsic calibrations 
(without dust); the third column provides the calibration for the 
reference model described above. The forth column provides 
the maximum value found in our models (the minimum value 
of calibrations corresponds to the case without dust). 

The calibrations are given assuming that the fraction of ion- 
izing photons absorbed by dust inside Hii regions (/d) is 0. 
Calibrations for different values /d can be obtained by dividing 
the value reported by 1 - fd- 

H/3. The first entry in the table is H/3 uncorrected for extinc- 
tion; this line may be the only hydrogen line detected in 
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the spectra of intermediate-redshift galaxies and, in that 
case there is no way to correct it for extinction. The cali- 
bration we provide here shows that for a typical SFR of a 
few Mg/yr, the inferred SFR may be underestimated by a 
factor of three. 

Hao. The Hao row refers to the Ha luminosity corrected for 
extinction using the Balmer decrement (ea. llOt . We notice 
that the decrement has been evaluated from comparisons of 
the observed ratio of the intensity of the emissions at Ha 
and HjS with the expected ratio. We did not attempt to sim- 
ulate a real measure of the Balmer decrement in a synthetic 
spectrum. Thus, our calibration assumes that one is able to 
correct the lines for the contribution of the underlying older 
populations. 

It is worth noticing that, due to the age selective extinction, 
the calibration may change by a factor of two, or even more 
because the extinction at Ha is underestimated, as demon- 
strated in Sect. 14.2.11 when the escape time is similar to the 
lifetime of ionizing stars. 

Another factor of uncertainty is provided by the orientation 
of the observed galaxy because the extinction from diffuse 
medium is characterized by a mixed geometry, which tends 
to produce a higher ratio between attenuation and redden- 
ing than for the screen geometry. 
Pa/?, Pa/3{). In Table |4] then Pa/3 follows, either uncoiTected or 
corrected (Pa/3o) for extinction, assuming that the intrinsic 
value of Pay6/Bry is 5.65 (and Tp^/rsiy = 2.95). In this 
case the extinction is lower and the SFR can be obtained 
with a large accuracy. However, it must be kept in mind 
that at these wavelengths uncertainties in the photometry 
of different wavebands may constitute the larger source of 
uncertainty in the estimated extinction (eg. Calzetti et al. 
[T996i l. 

Paa. The table continues with Paa uncorrected for extinction. 
This line sits at the short wavelength border of the K band 
(/l=18752A); like Hfi, it may be detected by ground obser- 
vations only in the spectra of intermediate redshift galaxies. 

Bry, Bi/3, Bra, Hua. The following entries in the table are 
Bry, BrjS, Bra and Hua uncorrected for extinction. As ex- 
pected, the extinction effects decrease when the wavelength 
increases, except for Hua (A = 12.369;um) that falls in the 
absorption feature of silicates. 

[On]3727. The table closes with the [Oii]3727 line uncor- 
rected for extinction (note we only use solar metallicity). 
The importance of this line is that in high redshift galax- 
ies this is the only bright line that, once redshifted, remains 
in the optical domain. Recently, Jansen et al. (1200 U inves- 
tigated the use of this line as SFR tracer, and found that 
SFR estimations by means of [O n]3727 agree with the Hao 
value, as long as the line intensity is corrected for extinc- 
tion and calibrated for metallicity. However, it is not pos- 
sible to perform these corrections without observing other 
lines that (if visible) are better estimators of SFR; therefore 
we conclude that [On] 3727 is a poor indicator of SFR. 



6. Infrared nebular metallic lines 

The IR spectral range contains numerous bright forbidden fine- 
structure lines that are excellent diagnostics of gas density, 
hardness of the exciting radiation field and abundance of im- 
portant elements. 



6.1. Density and ionizing spectrum hardness 
diagniostics 

Different IR fine structure transitions of the same ion have 
different sensitivity to collisional deexcitation. This can be 
used to identify the typical electron density of the emitting 
gas (Rubin et al. I1994> . Typical line pairs used for this pur- 
pose are [Arm]21.8jum/9//m, with a maximum sensitivity to 
the electron density at log(ne) = 4.7, [Ne iii]36/15.5;um (log(«e) 
= 4.7), [Om]52/88//m (log(ne) = 2.9), [Siii]19/33jum (log(«e) 
= 3.5), and [Nii]122/205//m (log(ne) = 1.8). Though colUsion 
strengths for these lines may need to be reexamined (Rubin 
et al. 200 Ij , it appears that they can trace the density of the 
emitting gas over a wide range, from hh = 10 cm"^ to ne - 
lO^cm"^, i.e. possibly encompassing the typical densities of 
star forming regions in normal galaxies as well as in compact 
obscured starbursts (Bressan et al. l2002> . 

The hardness of the ionizing spectrum and the average ion- 
ization parameter may be derived by analyzing ratios of line 
intensities in different ionization stages of suitable elements. 
Typical ratios are [Neiii]/[Nen], [Ar iii]/[Ar ii] and [Siv]/[Sm] 
(e.g. Giveon et al. l2002> . Our model can be used to obtain in- 
formation on the age of the ionizing stellar population and on 
the ionization parameter in single H ii regions or in starbursts of 
short duration. In other conditions, this method will only pro- 
vide the average values for these quantities. We will explore 
these topics in a forthcoming paper. 



6.2. Abundance determinations 

It is straightforward to extend to the infrared the method used 
for abundance determinations for optical lines. This method is 
based on the ratio between the sum of the intensities of the 
most abundant ions of an element, and that of a hydrogen re- 
combination line (typically the Bra line). It works as long as 
the current density is lower than the critical density for colli- 
sional deexcitation. 

In the near- and mid-infrared the elements to which we 
can apply the above method are neon ([Neiii]15.5/im and 
[Nen]12.8//m), argon ([Arm]9//m and [Arii]7//m) and sulphur 
([Siv]10.5yum, [Sra]18.7;um and [S ii]1.03//m). Suitable cali- 
brating relationships can be found in Giveon et al. J2002> and 
Verma et al. C2003j. 



In the far-infrared the method is applicable only to nitrogen 
([N ii] 122/im and [N iii]57jum). For these lines we have derived 
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Fig. 8. The relationship between nitrogen abundance (relative 
to hydrogen) and Nn (defined by ea. llSt for models with differ- 
ent densities and ionization parameters. The dashed line shows 
a linear fit to models with Z < 0.015. 



the following calibrations with respect to the Brc line (for dif- 
ferent assumptions on the electronic temperature Te): 



N 
H 



1.907 • 10-5 (L[N„]122 + 0.247L[Nl„]57) /iBr« 

forTe = 10000 K 

3.546 ■ 10-5 (L[N II] 122 + 0.231L[Niii]57) /i'Bra 

for Te = 5000 K 



(17) 



where N/H is the nitrogen abundance respect to hydrogen. The 
relevant transition data have been taken from Lennon & Burke 
l fT994j . Peng & Pradhan ( 1995) and Storey & Hummer ( 1995) . 

It is worth noticing that this method requires the measure- 
ment of the intrinsic intensity of a hydrogen recombination 
line. Thus, it is not applicable when these lines are not acces- 
sible, either because of instrumental limitation or for very high 
extinction. In this case, we suggest to adopt another indicator 
tightly correlated with the intrinsic intensity of hydrogen re- 
combination lines or, equivalently, with the SFR. The best SFR 
indicator we have found is the radio luminosity (Table |3}- To 
this purpose, we have run a set of models varying the metallic- 
ity of the gas (from Z=0.0008 to Z=0.05), the hydrogen den- 
sity, from 10 to 100 cm"^, and the ionization parameter The 
latter is a function of the mass of ionizing clusters M*, and 
of the filling factor e through eq. |8l We stress that the emis- 
sion lines resulting from quiescent star forming galaxies are 
the sum of the emission of H ii regions of different ages, so that 
there is not a single ionization parameter U at work. However, 
since all the ionization parameters scale linearly with the quan- 
tity m = (M*e^)'/^, we have changed m from 0.7 to 70 M^^ 
to simulate a similar range of variation of U. The fraction of 
ionizing photon absorbed by dust was assumed to be /d = 0.3. 

Fig. |8] shows the relationship in our models between the 
nitrogen abundance and the ratio N23 defined as 



We notice that the metallicity is well determined for values be- 
low solar For higher values the relation steepens considerably 
making the determination of the metal content quite uncertain. 

A linear fit to models in which the metallicity is lower than 
0.015 is plotted in Fig.|8j its analytical expression is 



log 



-18.21 + 1.203 logA^23 



(19) 



N23 - (i'[Nii]122 + 0.247L[Niii]57) /il.49GHz 



(18) 



where A^23 is given in hertz. 

We stress that this relation holds as long as the density is 
below the critical density for coUisional de-excitation of the 
[N 11] 122yL(m transition («e ~ 300 cm"-'). 

A similar calibration could be obtained by adopting the FIR 
luminosity instead of the radio luminosity. However, in the case 
of disk galaxies the relation between the FIR luminosity and the 
SFR is not as tight as that between the radio luminosity and the 
SFR (see Table|5J- This would introduce a significant scatter in 
the calibration. 

7. Discussion 

The first application of the model presented in this work is the 
analysis of attenuation in quiescent star-forming galaxies. 

These galaxies exhibit a poor correlation between the atten- 
uation in the UV band and Ha (B02, see also Fig.|5|i, with the 
former generally lower than what expected from the latter on 
the basis of the simple law Auv = l-6AH[r (eq.|5}. Our model 
explains this poor correlation in the context of age selective 
extinction. As the extinction in the UV is very high, the contri- 
bution to the observed UV flux from stars outside the molecular 
clouds is important even for relatively small optical thicknesses 
(tmc)- It is worth noticing that the fraction of young UV emit- 
ting stars outside MCs increases for decreasing escape time, 
while very young stars (which ionize the gas) spend most of 
their life inside MCs. Therefore, above a threshold value for 
Tmc which depends on the escape time, an increase of tmc pro- 
duces an attenuation in Ha larger than that in UV. 

Thus, neither observations nor modeling support the as- 
sumption of a constant relationship between the attenuation 
suffered by the continuum and the attenuation for the gas, at 
least for a current SFR smaller than 10 Mg/yr. As a conse- 
quence, the results of Calzetti (TWfl (A'^^^/A^^^ ~ 0.44) can- 
not be extrapolated from starbursts to disk star-forming galax- 
ies. 

Moreover, the extrapolation of the Ha attenuation to the 
UV through the previous reported simple law yields UV cor- 
rected fluxes larger than those expected from Ha corrected 
fluxes, as found by Sullivan et al. (2000 20011. These authors 
explain the result by introducing star formation histories that 
change rapidly with time and/or a more complex model of ex- 
tinction. On the other hand, our model explains this result as a 
natural consequence of the age selective extinction. 

An additional interesting issue is that the attenuation in 
Ha derived from the Balmer decrement can be underestimated, 
since in disk galaxies the escape time and the lifetime of very 
massive stars is often similar 

The comprehensive treatment of nebular and continuum 
emission allows to conclude that the assumption of a con- 
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Table 5. Comparison of some calibrations of SFR with results 
from previous works. Calibrations are referred to a Salpeter 
IMF between 0.1 and 100 M^. Reference: a: Kennicutt ( I1998> . 
b: Haarsma et al. (2000J. 



This work others 

UV 1.20 1.4 • 10-2' Moyr-'W-'Hz a 

IRg-iooo 4.63 4.5 ■ 10"" Moyr-'W-' a 

Ha 7.05 7.9 ■ 10^" Moyr-'W"' a 

1.49 GHz 79.6 123.3 ■ lO^^^ Moyr-'W-'Hz b 



slant extinction for all stellar populations is not satisfactory. 
However, the detailed models presented in this work indi- 
cate that Auv can be accurately estimated by using the ratio 
Ffir/F\j\/ (cfr. eas.ll3landll4>. This conclusion holds provided 
that the main contribution to FIR comes from the absorbed UV 
radiation of young stars, but it does not depend on the escape 
times and on the extinction properties of dust. The relation- 
ships we obtain are intermediate between those by Buat et al. 
jl999> and by Meurer et al. Cl999j . the latter derived for star- 
burst galaxies. 

By using our set of simulated disk galaxies, we calibrated 
different SFR estimators, collected in Tables |2]|3l and |3 Some 
calibrations may be compared with values obtained in the lit- 
erature. In particular we compare the results with Kennicutt 
i tTWHt (UV, IR and Ha) and Haarsma et al. (2000) (radio); 
these comparisons are summarized in Table|5j after a small cor- 
rection factor (1.16), which accounts for the slightly different 
IMF adopted, is applied to our estimates. 

There is a good agreement for Ha, UV and IR calibrations 
(note that Kennicutt I1998l uses the notation FIR for the total 
dust emission between 8 and 1000 /urn), while our calibration 
for the radio emission is a factor 1 .5 smaller than that proposed 
by Condon (1992) and Haarsma et al. (2000'), and a factor 1.3 
larger than the one quoted by Carilli (2002) . 

It is worth noticing that our calibrations give a ratio (eg. 1161) 
between the radio and FIR emissions, q - 2.32, very close to 
the observed value, and that we can explain in a natural way 
the scatter around the observed relation. In disk galaxies, this 
scatter is due to the variation of the fraction of stellar emission 
absorbed (and re-radiated) by dust. 

In the case of normal star forming galaxies, the model evi- 
dences that Ha, UV and even IR estimators of SFR are affected 
by important scatters. For what concerns Ha, this is due to the 
age selective extinction, that can produce an underestimate of 
the attenuation when it is derived by the B aimer decrement. 
Also IR emission by itself does not provide an accurate esti- 
mate of the SFR in disk galaxies, because the fraction of UV 
radiation not absorbed by the dust can be high. The combina- 
tion of UV and FIR luminosities provides a very good SFR 
estimator essentially because, even for a modest SFR, one re- 
covers in the FIR what is lost in the UV. 

Finally we have also discussed the utility of IR nebular 
lines as diagnostic tools for deriving the average properties of 
the environment such as number density, hardness of ionizing 
spectrum and gas metallicity. In particular we have provided 



a new calibration for the nitrogen abundance (N/H) as a func- 
tion of the intensities of the [Nii]122yum and [Niii]57jum, and 
Bra lines. When the latter line is missing or useless (e.g. when 
affected by strong extinction), we suggest the possible use of 
the radio luminosity as indicator of the ionizing flux Qn, that 
enters implicitely the above calibration. We provide a new cali- 
bration of the metallicity with the ratio N23 defined by equation 
[Tsl This new calibration will turn out to be particularly useful 
for the Herschel experiment. 



8. Summary 

Here we summarize the main points of this paper. 

1 . We have introduced the nebular emission calculations from 
CLOUDY into our spectrophotometric code GRASIL. The 
method which interfaces nebular emission computations 
with population synthesis is based on pre-computed li- 
braries of Hii region models. These libraries can be re- 
trieved from the GRASIL web site"* and easily used for 
other population synthesis codes. 

2. As a result, we can model the spectra of star-forming galax- 
ies from far-UV to radio wavelengths, including stellar ab- 
sorption features, nebular emission, dust and PAH emis- 
sion. 

3. We applied the model to study the attenuation in normal 
star forming galaxies. The poor correlation between the ex- 
tinction at UV and Ha found by B02 is interpreted in our 
model as a natural consequence of age selective extinction. 

4. Conversely, we found that Auv correlates quite tightly with 
the ratio Ffir/Fuv. 

5. We present new calibrations for the SFR using lines and 
continuum from UV to radio wavelengths, and study the 
reliability of each estimator. In particular, NIR H recombi- 
nation lines and radio luminosities are very accurate SFR 
estimators. 

6. Although UV-, Ha- (also corrected for extinction) and IR- 
derived SFR are plagued by significant uncertainties, the 
UV luminosity corrected by using the ratio Ffir/-Fuv is a 
robust SFR estimator 

7. We reproduce the observed value of the q ratio between ra- 
dio and FIR emission. Its scatter is ascribed to the variation 
of extinction between different objects. 

8. Finally, we discuss several methods and provide new cal- 
ibrations for the determination of metal abundances by 
means of infrared emission lines. 
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Appendix A: Analytical spectra 

Here we present the relationships used to describe the parame- 
ters of analytical spectra introduced in eq.^^s functions of Qh, 
Que, and Qq. 

Tiis is implicitly given by: 



1°° Fydv/hv By(Tn,)dy/hv Qn. 



(A.l) 



where VHe+ is the photoionization threshold frequency for He n. 
Once Tue is known, Aue is obtained from 



iHe 



Jvu, 



hv 



dv = Qhc 



(A.2) 



Furthermore, to compute Tu, we notice that there exists (for 
a wide range of age and metallicities) a correlation between Tu 
and QhJQh of the SSP specti-a, (Fig.lXH: 



3 ■ 104 + 4- 10^^ 



for ^ > 0.005 

Qh 2) 

4 ■ lO-* + 0.5 • lO-* ■ log(^) for ^ < 0.005 



We use this equation to derive Tu for given values of 
Qhs/Qh- Then Ah follows from 



(A.4) 



As for the values of Ani and of the temperature r,,!, we no- 
tice that the region of the spectrum with A > 912 A becomes 
important only for relatively low ionizing fluxes. We thus sim- 
ply obtained their values after minimization of the differences 
between the lines computed with full SSPs spectra and analyt- 
ical spectra: 



f Max(rH, 40000 K) for QhJQh > 0.005 



I 30000 K 

2.5Ah 
3. 5 Ah 



for QhJQh < 0.005 



for gHe/eH > 0.005 

for gHe/gH < 0.005 



(A.5) 



(A.6) 



In Fig. LA. 2! we compare some SSP spectra with the cor- 
responding analytical spectra. These comparisons suggest that 
analytical spectra can represent the Lyman continuum with rea- 
sonable precision. However, we stress that our goal is to get 
line emissions computed with analytical spectra which are as 
similar as possible to the ones computed with SSP. 

Appendix B: Average ionization parameter for a 
sphere 

The volume averaged ionization parameter for a shere is de- 
fined as: 



(U) 



-j 

Jo 



An p3 



(B.l) 



where is the Stromgren radius and U{r) is the ionization 
parameter of the gas at the distance r from the source. U{r) is 
defined by: 



Anr^riuc 



(B.2) 
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Fig. A.l. Temperature of ionizing continuum T^ as a function 
of Ghc/Qh for SSPs of difi'erent age and metallicity (points). 
The lines refer to the analytical fit. (ea. lA.3> . 
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Fig. A.2. UV and ionizing spectra for SSPs (Z=0.008) of dif- 
ferent ages (thick lines) compared with the corresponding ana- 
lytical spectra (thin lines) that have the same Q values than the 
SSP specti-a. 



where Q{r) is the number of ionizing photons per unit time that 
cross a spherical surphace at distance r. The variation of Qif) 
with the radius is given by: 



Integrating this equation we get: 
47r 

Q{r) = Qh - — r'4aB(H)e . 



(B.3) 



(B.4) 
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The Stromgren radius is defined as: 



(U) 



47r«2aB(H)« 

By substituting eas. lB.4l and lB.5l in eg. IB. II we get: 
I Qh rnHffB(H)e\ 3r^dr 

Jo 
3Gh 



(B.5) 
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